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ABSTRACT 

We investigate the evolution of Type Ib/c supernova (SN Ib/c) progenitors in close binary systems, 
using new evolutionary models that include the effects of rotation, with initial masses of 12 - 25 Mq 
for the primary components, and of single helium stars with initial masses of 2.8 - 20 M0. We 
find that, despite the impact of tidal interaction on the rotation of primary stars, the amount of 
angular momentum retained in the core at the presupernova stage in different binary model sequences 
converge to a value similar to those found in previous single star models. This amount is large enough 
to produce millisecond pulsars, but too small to produce magnetars or long gamma-ray bursts. We 
employ the most up-to-date estimate for the Wolf-Rayet mass loss rate, and its implications for SN 
Ib/c progenitors are discussed in detail. In terms of stellar structure, SN Ib/c progenitors in binary 
systems at solar metallicity are predicted to have a wide range of final masses up to about 7 Mq , with 
helium envelopes of A/hc — 0.16 — 1.5 M©. Our results indicate that, if the lack of helium lines in 
the spectra of SNe Ic were due to small amounts of helium (e.g. Mhc S 0.5), the distribution of both 
initial and final masses of SN Ic progenitors should be bimodal. Furthermore, we find that a thin 
hydrogen layer (0.001 Mq < Mh < 0.01 Mq) is expected to be present in many SN lb progenitors at 
the presupernova stage. We show that the presence of hydrogen, together with a rather thick helium 
envelope, can lead to a significant expansion of some SN Ib/c progenitors by the time of supernova 
explosion. This may have important consequences for the shock break-out and supernova light curve. 
We also argue that some SN progenitors with thin hydrogen layers produced via Case AB /B transfer 
might be related to Type lib supernova progenitors with relatively small radii of about 10 Rq. 
Subject headings: Stars:evolution, stars:rotation, stars:massive, binaries:close, supernovae: general 



1. INTRODUCTION 

It is generally believed that Type lb and Type Ic su- 
pernovae result from core collapse events of naked helium 
stars. The helium stars are thought to be produced by 
the loss of the hydrogen envelope, via stellar winds mass 
loss from massive single stars or via mass transfers in 
close binary systems. 

According to recent stellar models adopting the 
most up-to-date st ellar winds i nass loss rates 
jMevnet & Maeder ' '2003', '2 0051: lEldridge fc Vinki 
[2006; Limongi & Chieffi 2006:~ IGeorgv et al. I 12009^ 
the final masses (Mf) of helium stars produced by 
mass-losing single stars appear to be too massive 
to produce typical SNe Ib/c (i.e., Mf > 10 Mq at 
solar metallicity). Although the limiting mass for BH 
formation is not yet well determined, given their high 
binding energy, such massive progenitors of A/f > 10 Mq 
are likely to form black holes (BHs), p roducing fain t 
supernovae or no supernova at ah (cf. iFrverl Il999h . 
Although very bright SNe Ib/c like SN 1998bw could 
be produced from such massive helium stars if, for ex- 
ample, pow ered by rapid rotation fe.g.. I Wooslev I [19931 : 
iBurrowset al. 2007[). such events are shown to be rare 
(e g.. iPodsiadlowski et 3171120041: iGuetta fc Delia Valle I 
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By contrast, helium stars with a wide range of masses 
(2.0 Mq < A/hc < 25 Mq) can be made from 12. ..60 Mq 
primary components in close binary systems via the so- 
called Case A/B mass transfer[3 Many of them may end 
their life as bright SNe Ib/c leaving neutron stars as rem- 
nants, if their final masses are less than about 7-10 Mq. 
Population studies indeed show that close binary stars 
can produce a sufficient number of SNe Ibc to explain 
their observed rate, without the need of invoking single 
star progeni tors (e.g. , IPodsia dlowski. Joss fc Hsul 

19921: Ide Donder & VanbcverenJ 119981 : 

Eldridge. Izzard k Touti i2008il Therefore, it is 
most likely that the majority of typical SNe Ibc are 
produced in binary systems. 

The observational evidence for the connection be- 
tween SNe Ibc and long gamma-ray bursts (GRBs) 
has particularly motivated many observational studies 
to b etter understand SNe Ib c since the last decade 
(see iWooslev fc BloonTI 120061 for a review). The- 
oretical stellar models of SNe Ibc progenitors are 
thus highly required nowadays. The most compre- 
hensive studies on the detailed characteristics of SNe 

Case A, B or C mass transfer denotes mass transfer from the 
primary star during core hydrogen burning, helium core contrac- 
tion/beginning of core helium burning, or core helium burning and 
later stages, respectively. On the other hand, if mass transfer oc- 
curs during helium core contraction/beginning of core helium burn- 
ing from a star that has already undergone Case A mass transfer, 
such a mass transfer phase is called Case AB. Case ABB or Case 
BB mass transfer denotes mass transfer from the primary star dur- 
ing core helium burning and/or later stages, which has already 
undergone Case AB or Case B mass transfer, respectively. 
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Ibc progenitors in binary system s were conducted by 
iWooslev. Langer fc Weaverl (flOOSl) (hereafter, WLW95) 
using mass-losing pure h elium star models, and by 
iWellstein fc Langer I (jl999( ) (hereafter, WL99) using self- 
consistent binary star models. Although more recent the- 
oretical studies on SNe Ibc progenitors in binary systems 
can be found in the literature, they have been focused 
on long GRB progenitors or stellar populations, rather 
than on the detailed nature of typical SNe Ibc progenitors 
(e.g., 'Brown e t al. ll20p( J: 'Izzard. Ramirez-Ruiz & Tout ' 
2005a; Cant iello et al. 2007: 
' |200^ fPetmers et al. I I200I 

\2om . 
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In this paper, we revisit the problem of SNe Ibc progen- 
itors in close binary systems using both binary star and 
single helium star models up to the neon burning stage, 
with updated physics of two important ingredients. One 
is rotation, which was not considered in WLW95 and 
WL99, and the other is the mass loss rate of Wolf-Rayet 
(WR) stars (Sect. [2]). 

This paper is organized as follows. In Sect.[2l we briefly 
review recent developments of stellar evolution models 
regarding the effects of rotation and the WR star mass 
loss rate, arguing for the need of updated physics in bi- 
nary star models. Our adopted physical assumptions and 
numerical method are discussed in Sect. [3] In the fol- 
lowing section (Sec. |4]), using our binary star evolution 
models including the effect of rotation and the transport 
of angular momentum due to hydrodynamic instabilities 
and magnetic torques, we explore the role of tidal inter- 
action and mass transfer in the redistribution of angu- 
lar momentum in primary stars. In Sec. [Sj the nature 
of SNe Ibc progenitors is investigated in terms of final 
masses, masses of helium and hydrogen layers, radii and 
mass loss rates at the presupernova stage, assuming these 
properties do not significantly change from neon burning 
to core collapse. For this purpose, we also present mass- 
losing single helium star models as a complement to our 
binary star models, given that the parameter space ex- 
plored with our binary model sequences is limited. We 
conclude the paper by discussing observational implica- 
tions of our results, in Sect. El 

2. ROTATION AND WOLF-RAYET WINDS 
2.1. Rotation 

Rotation has particular roles in the evolution 
massive stars as it changes the stellar structure, 
induces chemical mixing and enhances ma ss loss 
due to stellar winds (iMaeder fc Mevnet I 120001: 
iHeger. Langer fc Wooslevl I2000D . During the last 
decade, several authors have calculated massive star 
models up to the pre-supernova stage, considering 
the redistribution of angular momentum and chemical 
species due to rotationally induced hydrodynamic 
instabilities, such as Edd ington-Sweet circulations, 
and t he shear instability (iHeger Langer fc Wooslev I 



[Mot iHirschi. Mevnet fc Maeder I l2004f ). Afthough 
these models could explain some observational aspects 
such as surface abundances of CNO elements of mas- 
sive stars and Wolf-Rayet ( WR) star p opulatio ns at 
different metallicities (e.g.. | Maeder fc "M evnet 120001 : 
IHeger fc Langer 1 120001: iMevnet fc Maeder 1 i2005i) . their 
adopted angular momentum transport mechanisms 



turned out to be too inefficient to explain the ob served 
spin rates of stellar remnants ()Suiis et al. 1 120081 ). I.e., 
their models predict nearly two orders of magnitude 
higher spin rates of white dwarfs and young neutron 
stars than the observed ones. These models also imply 
that almost all WR stars can retain enough angular 
momentum to produce GRBs, either by magnetar 
or collapsar formation depending on the final mass. 
This gives a nearly 1000 times higher ratio of GRBs 
to SNe than the observationall y imp lied value. On 
the other hand, ISpruit ' (199^ I2002D suggested that 
magnetic torques resulting from dynamo actions in 
differentially rotating radiative layers (the so-called 
Spruit-Tayler dynamo) should be the dominant angular 
momentum transport mechanism compared to the pure 
hydrodynamic instabilities. Recent magnetic models 

that adopt the Spr uit-Tayl er dynamo according to 

the prescription by iSpruit I (|2002i ) are indeed more 
consistent with observ ations in terms of the spin rates 
of stellar remna nts (jHeger. Wooslev fc Spruit I 120051 : 
iSuiis et al. II2008D . Magnetic models also better explain 
the fact that long GRBs are rare events compared 
to normal core collapse supernovae. Although the 

Spruit-Tayler dyna mo mechanism is still subje ct to 

many uncertainties (iDenissenkov fc Pinsonneaultl 120071 : 
iZahn. Brun fc Mathis I I2007D . these recent studies 
indicate that an efficient angular momentum transport 
mechanism comparable to what the Spruit-Tayler dy- 
namo predicts is needed to understand the observations. 

The above discussion is based on single star models, 
and the role of rotation in the evolution of massive bi- 
nary stars remained relatively unexplored. Massive stars 
in close binary systems are supposed to experience an ex- 
change of mass and angular momentum via mass trans- 
fer and tidal interaction, and thus the evolutio n of binary 
stars i s more complex than that of single star s . IWellstein I 
(|2001D and iLanger. Wellstein fc Petrovic I (|2003D pre- 
sented, for the first time, self-consistent calculations of 
massive binary star evolution models including many rel- 
evant effects of rotation and binary interactions: the 
change of the stellar structure due to the centrifu- 
gal force, transport of angular momentum and chem- 
ical species due to rotationally induced hydrodynamic 
instabilities, tidal interaction, transfer of mass from 
the primary, accretion of mass and angular momentum 
of the secondary, and resulting changes of the orbit. 
Their work indicates that the secondary can be eas- 
ily spun up by mass accretion even up to critical rota- 
tion, thus modulating the mass accretion efficiency by 
the interplay of the enhanced mass loss due to rota- 
tion from the spun-up secondary and the mass trans- 
fer from the primary. This effect provides important 
clues to better understand the evolutionary paths of 
some observe d X-ray and WR star binary systems as 
discussed by iLanger. Wellstein fc Petrovic I ([2003) and 
iPetrovic. Langer fc van der Hucht I (l2005bD . Their non- 
magnetic models also show that massive stars may end 
up with different core spin rates depending on the his- 
tory of mass loss/gain during binary evolution, that could 
be related to the observatio nal diversity of c ore-col lapse 
supernovae. Mo r e rece ntly, IPetrovic et al. I ()2005al ) and 
iCantiello et al. I (|2007D included the Spruit-Tayler dy- 
namo in their binary models, and discussed possible 
evolutionary paths of massive binary stars towards long 
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Fig. 1. — Wolf-Rayet mass loss rates from helium stars on the 
zero-age main sequence at s olar meta llicity. The mass loss rates of 
Hamann et al. (see Eq. [!}, ILan ger 1 fl989) and Nugis & Lamers I 
([2000) are given by dotted, dashed and solid lines, respectively. 
The dotted-dashed and three dotted-dashed lines denote the 
Hamann et al. rates divided by a factor of 5 and 10, respectively 
(i.e., /wH, = 5 and 10 in Eq. [T]). 

GRBs. In this paper, we present new evolutionary calcu- 
lations of magnetic (i.e., the Spruit- Tayler dynamo is in- 
cluded) massive binary stars for a large parameter space, 
focusing on the evolution of the primary stars to investi- 
gate the nature of typical SNe Ibc progenitors. 

2.2. Mass loss due to Wolf-Rayet winds 



WLW95 employed the 
loss rate of WR stars 
their helium star models, 
of iHamann . Schonber ner 



m ass-depe nd ent mass 
of ILanger I (pSOh for 
and WL99 used those 
and 



IHaman n. Kocs terke fc W cssolowski 
studies concluded that the final 



Both 

masses of SNe Ibc 
progenitors in binary systems should converge to a 
limited range of 2.2 - 3.6 Mq, even for an initial helium 
star mass of 20 Mq. Later developments of the WR 
wind theory considering dumpies pointed out, however, 
that the WR mass loss rates used by WLW95 and WL99 
are si gnificantly overes timated (jHamann fc Koesterke I 
[1991 ikugis fc La^^Ierel [2OOOI ). For example. Fig. [D 
shows that, for helium stars on th e zero age main 
sequen ce, the mass loss rate given by iNugis fc Lamersi 
(|200C0 is almost an order-of-magnitude lower than that 
of Hamann et al. used in WL99. The most recent 
theoretical models of WR winds by Vink fc de Koterl 
(jlOOS) and Gra fener fc Hamann I (j^08) also give WR 
mass loss rates compatible to the Nu gis & Lamers rate 
(cf. Figure 1 in lYoon fc Langer I [2001 . 

The reduced WR mass loss rates have been consid- 
ered in m any recent massive single star models (e.g., 
iMevnet fc Macdcr 2003, 2005; Eldridgc & Vink 2006), 
suggesting large final masses of WR stars at the presuper- 
nova stage (Mf > 1 at solar metallicity) compare d 
to - 4 Mq found bv lWooslev. Langer fc Weaver I (jl993[ ). 
As argued in the introduction, this leads to the interest- 
ing conclusion that most core collapse events occurring in 
single WR stars should not produce bright supernovae, 
rendering the binary star channel even more important 
for SNe Ibc. Surprisingly, the effect of the new WR 
mass loss rate on SNe Ibc progenitors in close binary sys- 
tems has not been much discussed or la rgely overlooked. 
To our knowledge. [Pols fc Dewi[ (|200^ are the only au- 
thors who addressed this issue in detail. Although WL99 



also presented some model sequences with a reduced WR 
mass loss rate, their adopted reduction factor was rather 
modest (2 times smaller than the Hamann et al. rate). 
Here we consider the Hamman et al. WR mass loss rate 
reduced by factors of 5 and 10, which reflects the most 
recent result as discussed above. 

3. NUMERICAL METHOD AND PHYSICAL ASSUMPTIONS 

We have used the sam e stellar evolution code as 
in ICantiello et al. I (|2007[ ). that follows the simultane- 
ous evolution of the two stellar components of a bi- 
nary system. The effect of the centrifug al force on 
stellar structure is considered following lEndal fc Sofia I 
(|1976l ). The transport of angular momentum and 
chemical elements is treated as diffusion consider- 
ing Eddington Sweet circulations, the shear instabil- 
ity, the Goldreich-Schubert-Fricke inst ability and the 
Spruit- Tayler dynamo, as described in iPetrovic et al. I 
(|2005aD . We use g .sRMT = 1.0 for the semi-co nvection 
parameter (ILanger. Fricke fc Sugimoto I 119830 . as in 
lYoon. Langer fc Norman I ( 2006f ). 

The stellar winds mass loss rate is determined fol- 
lowing (jKudritzki et al. 1I1989D with a metallicity scal- 
ing of (Z/Zq)0-69 (|Vink. de Koter fc Lamersi [200l for 
the main sequence phase. WR wind mass loss rates are 
computed according to Hamman et al. with a correc- 
tion factor /wR, and with a metallicity dependence of 
M oc (Zi„it/Z0)O-86 (|Vink fc de Koter ll2005l) : 



log 



WR 



Mq yr- 



= -11.95-1- 1.5 log L/Lq - 2.85Xs 

+0.86 log(Zi„it/Zo) - log /wR 
for logL/L© > 4.5, (1) 
= -35.8-l-6.81ogL/L0-2.85Xs 

+O.861og(Zi„it/Z0)-log/wR 
for logL/L0 < 4.5 . 

We use /wR = 5 or 10 in most model sequences, which 
means the WR mass loss rate by Hamman et al. is low- 
ered 5 or 10 times, to consider the most recent estimates 
(see Fig. [1]). The enhancement of stellar winds mass loss 
due to the centrifugal force is considered by 



M 



M{Vrot = 0) 



1 



1 - n 



0.43 



,0.5- 



M 



(2) 



where fl — Wrot / ^crit wit h fcrit — \/GM (l — T)/R and F 
is the Eddington factor (ILanger II1998D . Here, to prevent 
a singularity that may occur as Wrot approaches Ucrit, the 
mass loss rate is limited to M < 0.5M/rKH, where tkh 
is the thermal time scale of the star. 

The binary orbit is assumed to be circular, and the 
Roche lobe ra dius is det e rmine d according to the ap- 
proximation of lEggleton I (|1983[ ). The mass loss rate of 
the Roche-lobe filling component through the first La- 
grangian point is implicitly computed using the method 
given by Rittcr (1988). The equation of motion of 
a test particle is numerically solved to calculate the 
amount of angular momentum of the accreted matter 
if the transfered matter directly hits the seco ndary star, 
and t he Keplerian value is assumed otherwise (iWellstein I 
l2001f) . The change of the orbital period due to mass 
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transfe r an d stellar wind mass loss is cons idered accord- 
ine; to [ Pod siadlo wski. Joss fc Hsul ()1992l ). We follow 
iBrooksbiaw fc Tavani I ()1993[ ) to determine the amount 
of the specific angular momentum carried away from the 
orbit by stellar winds. 

Tidal sy nchro nization is considered fo l lowin g; 
iWellstein I (|200l (see also iDetmers et aFI l2008l ). 
We assume a synchro nization time scale according to 
iTassoul I (|1987l I2000D who considered tidally driven 
meridional circulations as the main mechanism for tidal 
dissipation: 



1/4 



1.44 X 10^ '' /Lq 

9/8 /^x 33/8 



(3) 



where q denotes the mass ratio and d the orbital sep- 
aration. This prescripti on gives a m uch shorter time 
scale than that given by iZahn I (|1977| ) . Given that the 
physics of tidal dissip ation is much debated in the liter- 

to 



ature (jLanger 



L oissip; 

Mm, 



we introduce a parameter /g, 



investigate how an extremely fast/slow synchronization 
may influence the results. In most cases, however, we use 

/sync = 1- ^ 

A fe w sequences are also computed with Tgync of IZahn I 
(flQTl for comparison: 



' sync 



fGM\ 



1/2 



MR^ [R 



17/2 



(4) 



where / is the moment of the star, and E2 a constant 
measuring the coupling between the tidal potenti al and 
the gravity mode. Using the data of Table 1 in IZahn I 
we constructed a fitting formula for as the 
following: 



= 10- 



-1.37 



Rc 



R 



(5) 



where -Rconv is the radius of the convective core. Note 
that both prescriptions by Tassoul and Zahn are not ap- 
propriate for a star with a convective envelope 0. How- 
ever, the role of tidal synchronization is significant only 
on the main sequence, and not important in late evolu- 
tionary stages as discussed below. 

We computed 45 model sequences for initial masses of 
the primary star mostly from 12 to 25 M0 at two differ- 
ent metallicities {Z = 0.02 and 0.004), for different mass 
ratios, initial orbital periods, and WR mass loss rates, as 
summarized in Table. [TJ The initial rotational velocity 
at the equatorial surface of each star is set to be 20% 
of the Keplerivan value. We could not calculate more 
massive systems because of a numerical difficulty encoun- 
tered during the mass transfer phases, except for Seq. 26 
where a primary star of 60 Mq is considered with a rather 
large WR mass loss rate (i.e., /wr = 3). The adopted 
initial orbital periods corresponds either to Case A or to 



^ O n the other hand, note that a recent study bv lToledano et al. I 
1)20071 ) suggests that intermediate mass main sequence stars follow 
the Zahn's synchronization time scale for convective stars. 
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Fig. 2. — Evolution of the internal structure of the primary star 



in Seq. 14 (Mi^i^it = 18 MQ.Mj.init = 17 Mq and P;,, 



4 day) 



from ZAMS to the neon burning phase. The hatched lines and 
the red dots denote convective layers and semi-convective layers, 
respectively. The different shades give the nuclear energy genera- 
tion rate, for which the scale is shown on the right hand side. The 
surface of the star is marked by the topmost solid line. 

Case B mass transfer. In the present study, we do not 
consider Case C systems, but briefly discuss the possible 
outcomes of Case C mass transfer in Sect. 16.41 The evo- 
lution of the primary stars is followed up to neon burning 
in most cases. 

We also present non-rotating single helium star models 
to discuss SNe Ibc progenitors in binary systems with 
initial masses larger than 25 M©, and also to compare 
them with binary star models (Sect. [S]). 

4. REDISTRIBUTION OF ANGULAR MOMENTUM IN 
PRIMARY STARS 

In this section, we focus our discussion on the evo- 
lution of primary stars and investigate whether binary 
evolution via Case A or Case B mass transfer could lead 
to diverse pre-coUapse conditions of SNe Ibc in terms 
of the amount of core angular momentum. Although the 
evolution of mass- accreting secondary stars is a matter of 
extreme interest as disc usse d inlBraun fc Langer I (119951) , 
iPetrovic et aFI ()2005aD and lCantiello et al. I ()2007D . it is 
beyond the scope of this paper. 

Here, we first present some results including the Spruit- 
Tayler dynamo with our fiducial assumption on synchro- 
nization time (i.e., /sync = 1)) showing that the final 
amount of angular momentum in the core of the primary 
star is not much affected by different histories of mass 
loss (i.e.. Case AB or Case B; Sect. 14.1] ). Then, we dis- 
cuss the influences of different assumptions on tidal syn- 
chronization and transport process of angular momen- 
tum (Sect. 1121). 

4.1. Fiducial Models 

4.1.1. Evolution with Case A and AB mass transfers 

The evolution of the primary star in a close binary 
system is characterized by the rapid loss of mass due to 
Roche-lobe overflow. As an example, the evolution of 
the primary star in Seq. 14 is described in Figs. [Hand |31 
where our fiducial value of /sync = 1 is adopted, includ- 
ing the Spruit-Tayler dyanamo. The binary system ini- 
tially consists of a 18 star and a 17 Mq star in a 
4 day orbit. Mass transfer starts at t — 8.09 x 10^ yr, 
when the helium mass fraction in the hydrogen burn- 
ing core has increased to 0.94. The mass transfer rate 
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Fig. 3. — Upper panel: Evolutionary track of the primary star 
in Seq. 14 (Mi^i^it = 18 Ma.Mj.init = 17 Mq and Pi^n = 4 day) 
in the HR diagram. The filled circles on the track mark different 
evolutionary epochs as the following. 1: ZAMS, 2: beginning of the 
Case A mass transfer, 3: end of the Case A mass transfer, 4: core 
hydrogen exhaustion, 5: beginning of the Case AB mass transfer, 
6: end of the Case AB mass transfer, 7: core helium exhaustion, 
8: neon burning (end of calculation). Lower panel: Mass transfer 
rates from the primary star (solid line) and mass accretion rates 
onto the secondary star (dashed line) during the Case A and AB 
transfers as a function of the primary star mass. 

rises up to 8 x lO"** MQyr~^, which roughly corresponds 
to Mi/tkh,i where Mi and tkh,i denote the mass and 
the Kelvin-Hehiiohz time scale of the primary star, re- 
spectively. The primary mass decreases to 7.5 Mq by 
the end of the Case A transfer (see Fig. U). The sec- 
ond Roche-lobe overflow begins at t = 8.513 x 10^ yr 
when the envelope of the primary star expands due to 
hydrogen shell burning during the helium core contrac- 
tion phase (Case AB mass transfer). The primary star 
loses most of the hydrogen envelope as a result, expos- 
ing its helium core of 3.95 Mq having a small amount of 
hydrogen (Mh — 0.04 Mq) in the outermost layers, as 
shown in the third panel of Fig. |4l 

Although the star remains compact (i? < 0.9 Rq) dur- 
ing core helium burning, helium shell burning activated 
after core helium exhaustion leads to the expansion of the 
envelope up to ~ 12 Rq (see Fig. [3]) during core carbon 
burning. A Case ABB mass transfer does not occur, how- 
ever, due to the large orbital separation {A 121 Rq) 
at this stage, while it does occur in many other sequences. 
The final mass at the end of the calculation (neon burn- 
ing) is 3.79 Mq. The mass of hydrogen decreases to 
0.0015 Mq at the end, and the remaning mass of helium 
is 1.49 Mq, as shown in the last panel of Fig. S) The star 
is likely to eventually explode as a Type lb supernova 
given the rather thick helium envelope with a very thin 
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Fig. 4. — Chemical composition of the primary star in Seq. 14 as 
a function of the mass coordinate, at different evolutionary epochs. 
First Panel: Core H burning (right before the Case A mass transfer 
phase) Second Panel: Core H burning (right after the Case A mass 
transfer phase) Third Panel: Helium burning (right after the Case 
AB mass transfer phase) Last Panel: Neon burning 
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hydrogen layer, but it might also appear as Type lib if 
the supernova were found within several days after the 
explosion (see Sect. 15. 3p . 

The lower panel of Fig. [3] shows that the mass trans- 
fer is not conservative. When the secondary star reaches 
critical rotation as a result of the accretion of angular 
momentum, the stellar wind mass loss rate increases 
so drasti cally as to prevent efficient mass accumula- 
tion (see iPetrovic. Langer fc van der Hucht I [2005R for 
a more detailed discussion on this effect). The ratio of 
the accreted mass in the secondary star to the trans- 
ferred mass from the primary star is about 0.83 during 
the Case A mass transfer, and 0.41 during the Case B 
mass transfer. The previou s calcu lations by WL99 and 
iWellstein. Langer &: Brauni ()2001[ ) show that the non- 
conservative mass transfer leads to a shorter orbit than 
in the case of conservative mass transfer, in general. This 
effect should be kept in mind in the following discussion 
on the evolution of the binary orbit and its consequences. 

The evolution of rotation in the primary star is shown 
in Figs. [5l |6] and [71 As shown in Fig. [5l the total an- 
gular momentum of the primary star rapidly decreases 
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in the beginning as a result of tidal interaction, until 
the star is completely synchronized with the orbit when 
t ~ 2 X 10^ yr. Fig. [5] shows the evolution of the an- 
gular velocity inside the primary star from ZAMS un- 
til it reaches complete synchronization. Note that the 
star keeps rotating almost rigidly, even though the tidal 
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interaction causes redistribution of anguiar momentum 
from the outermost l ayers of tlie star |Wellstcin 2001; 
iDetmers et al. I [20081 ) . This results from the very short 
time scale of the transport of angular momentum in the 
star (about 300 years) due mainly to convection and the 
Spruit-Tayler dynamo, in the convective and radiative 
layers respectively. 

Later on, the angular velocity of the primary star still 
remains coupled with the orbital motion, until decou- 
pling starts during Case B mass transfer (see below) . In- 
terestingly the total angular momentum of the primary 
star gradually increases from the initial synchronization 
until the onset of Case A mass transfer (Fig. [5]). This is 
because the star significantly expands, while the change 



of the orbital separation remains small durning this pe- 
riod, as shown in the third and last panels in Fig. [T) 
However, it rapidly decreases again during the Case A 
and AB mass transfer phases, as explained below. 

Fig.[7]shows that angular momentum in the core of the 
primary star is mostly removed during the mass transfer 
phases (Case A and Case AB). It should be noted, how- 
ever, that the mechanism for braking the core is different 
for each case. During the Case A mass transfer, the syn- 
chronization time scale according to Eq. ([3]) remains very 
short (10'^ yr) compared to the mass transfer time scale 
( 10^ yr). The decrease of the core angular momentum 
during the Case A mass transfer results from the syn- 
chronization that occurs even when the orbit is rapidly 
widened due to mass exchange (Fig. [7]). During the Case 
AB mass transfer, the further increase of the orbital sep- 
aration significantly weakens the role of synchronization 
for the redistribution of angular momentum. However, 
both the Spruit-Tayler dynamo and mass loss lead to 
rapid braking of the thermally contracting helium core. 
Further significant core-braking by the Spruit-Tayler dy- 
namo occurs during the CO core contraction phase, and 
the mean specific angular momentum of the innermost 
1.4 M0 becomes about 2.5 x 10^"* cni^ s^^ at the neon 
burning phase. 

4.1.2. Evolution with Case B mass transfer 

In Seq. 20, the initial masses of the stellar compo- 
nents are the same as in Seq. 14 but the initial period is 
large enough that the first mass transfer occurs during 
the helium core contraction phase (Case B mass trans- 
fer). Avoiding Case A mass transfer, the primary star 
retains more angular momentum in the core at the end 
of main sequence than in Seq. 14. However, the core loses 
more angular momentum during the helium core contrac- 
tion phase than in Seq. 14, as stronger magnetic torques 
are exerted, mainly due to the more massive envelope 
(Fig. H]). At neon burning, < J1.4 >= 2.6 x 10^^ cm^ 
is obtained, which is similar to that in Seq. 14. 

4.2. Non-fiducial Models 
4.2.1. Influence of the synchronization time scale 

The effect of synchronization is negligible in Seq. 15, 
where /sync = 10^ is adopted. As shown in Fig. |9l a 
rather rapid decrease of < jcoro > occurs during and 
after the Case A mass transfer, which is a combined ef- 
fect of mass loss and magnetic torques: mass loss car- 
ries away angular momentum from the envelope, and 
magnetic torques brake the core rotation subsequently. 
The core is further slowed down during the helium core 
and CO core contraction phases, resulting in < ji.4 >= 
2.6 X 10^^ cm^ s~^ at the neon burning phase. Note that 
this value is very close to that in Seq. 14. 

We find that the result with Zahn's prescription for 
synchronization (Seq. 17) is not much different from that 
of Seq. 15 where synchronization is negligible. The syn- 
chronization time scale according to Zahn is sensitive to 
the ratio of the convective core size to the stellar radius 



sync 



(i?co„v/i^star)^ scc Eqs. ^ and dSl)). This ra- 
tio continuously decreases as the star evolves, and thus 
T'sync contiuuously increases to such an extent that the ef- 
fect of synchronization can be ignored when Case A mass 
transfer starts. When /sync = 0.01 with Tassoul's pre- 
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scription is used (Seq. 16), on the other hand, synchro- 
nization becomes important even after the Case AB mass 
transfer phase, and the primary is further spun down by 
tidal interaction, giving < ji.4 >= 6 10^^ x cm^ s^^ at 
the neon burning phase. 

4.2.2. Non-magnetic model 

In Seq. 18 where the Spruit-Tayler dynamo is not in- 
cluded, the core is spun down due to synchronization 
during the Case A mass transfer phase (Fig. |9]) as in 
the corresponding magnetic case (Seq. 14; Fig. [71. The 
spin-orbit couphng becomes significantly weakened as the 
orbit widens after the Case A mass transfer phase. De- 
spite a significant amount of mass is lost via Case AB 
mass transfer, the core in the primary star retains most 
of the remaining angular momentum in the following evo- 
lutionary stages. This is because the chemical gradi- 
ent across the boundary between the helium core and 
the hydrogen envelope effecti vely prohibits the trans- 
port of angular momentuni (cf.lMevnet fc Maeder II1997I : 
iHeger. Langer fc Wooslev 1120001 ). The core angular mo- 
mentum at neon burning is thus about 10 times larger 
(< J1.4 >= 3.57 10^^ cm^ s~^) than in the corresponding 
magnetic case. 

4.3. Discussion 

As shown in the above examples, in the model se- 
quences with the Spruit-Tayler dynamo, all of the pri- 
mary stars retain similar amounts of angular momen- 
tum (a few 10^** cm^ s~^) in the innermost 1.4 M© at 
neon burning regardless of the detailed history of mass 
transfer, unless synchronization is extremely fast as in 
Seq. 16 (see Table [T]). According to the Spruit-Tayler 
dynamo, magnetic torques exerted to the core become 
stronger with a higher spin rate, a larger degree of differ- 
ential rotation between the core and the envelope, and a 
heavier radiative envelope. Therefore, although winds 
or Roche-lobe overflows reduce the size of the hydro- 
gen envelope and remove angular momentum from the 
star, this in turn weakens the torque exerted to the core, 
and vice versa. The remarkable convergence of < J1.4 > 
to a few lO^** cm^ s~^ in our model sequences, even for 
different wind parameters and metallicities as shown in 
Table [U can be explained by this self-regulating nature 
of the Spruit-Tayler dynamo. This result indicates that 
not much diversity is expected in SNe Ibc progenitors 
produced via Case A or Case AB/B mass transfer, in 
terms of rotation: most of SNe Ibc of a similar progen- 
itor mass may leave neutron stars with a similar spin 
rate. However, other types of binary interactions still 
may lead to various final rotation periods in SN pro- 
genitors (sec Brow n et aTl 120001 : iCantiello et al. 1 120071 : 
iPodsia dlowski ct aLl l201(]l for such examples). 

5. THE NATURE OF SN IBC PROGENITORS 

In the literature, the detailed history of binary interac- 
tions such as mass transfer and tidal interaction is often 
neglected, and only the evolution of pure helium stars 
to discuss SNe Ibc prog enitors in binary sy stems is in- 
vestigated (e.g. WLW95: lPols fc Dewi |[2002l ). Although 
this approach is very useful, our binary star models in- 
dicate that complications resulting from binary interac- 
tions have important consequences in SN Ibc progenitors. 
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circles denote the prediction from our binary star models. The 
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To address this issue, here we also present evolutionary 
models of non-rotating single helium stars with initial 
masses of 2.8 - 20 Mq at Z = 0.02 for comparison with 
binary star models, as presented in Table[2j This also en- 
ables us to study SN Ibc progenitors with initial masses 
higher than 25 Mq, which are lacking in our binary model 
sequences. These helium star models were calculated up 
to neon burning, with the WR mass loss rate given in 
Eq. ([T]) with /wR = 5 or 10, as in the binary models. 
In the following, we focus our discussion mostly on the 
results with /wr = 5, while the influence of the wind 
parameter is briefly discussed. 

5.1. Final mass 

The final mass of the primary stars in a close binary 
systems is largely determined by the winds mass loss dur- 
ing core helium burning, if the initial mass (Mhc.i; i.e., 
the mass right after Case AB or Case B mass transfer) 
is significantly larger than about 3.0 Mq. However, for 
a less massive helium star, the expansion of the helium 
envelope becomes so dramatic during CO core contrac- 
tion, and/or during core carbon burning, that it can lead 
to another mass transfer phase: Case ABB or Case BB 
(e.g., compare the final radius of the 2.8 Mq helium star 
model with those of more massive helium star models in 
Tabic El). 

The impact of this mass transfer phase becomes more 
important for a less massive helium star. In Seqs. 1, 2 
and 4 where MhcI — 2.1 - 2.3 Mq, the rapid loss of mass 
during Case ABB/BB transfer reduces the total mass 
of the primary stars to such an extent that they may 
not explode as supernovae, but die as white dwarfs. In 
Seq. 3, where Mho,! — 2.7 Mq, the carbon-oxygen core 
can grow beyond the Chandrasekhar limit despite the 
significant loss of mass (i.e., about 1 Mq) during the Case 
BB phase. The remaining helium mass in the envelop is 
only about 0.18 Mq at core carbon exhaustion. For the 
other sequences at solar metalhcity (Seqs. 5, 8, 9, 10, 11, 
13, 21) where MhcI > 3.0 Mq, the amounts of mass loss 
via Case ABB/BB transfer until core neon burning are 
rather moderate, varying from 0.02 to 0.1 Mq. 

The WR wind mass loss rate has significant conse- 
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Fig. 11. — The amount of hydrogen in SN Ibc progenitor models 
as a function of the final mass. The filled triangles and circles 
denote the results of our binary star models at Z =7,q with /wR = 
5 and 10, respectively. The fille squares give the results with the 
SMC metallicity models. 

quences in systems where Afnc.i > 3.0 M©. The pre- 
vious study by WLW95 concluded that the final masses 
of single helium stars with initial masses of 3 - 20 M© 
should converge to ~ 2.2 - 3.5 M©. As expected from the 
lower WR mass loss rate adopted in this study, our he- 
lium star models give a significantly wider range of final 
masses (Mf) at the given metallicity of 0.02: M{ = 2.9 - 
7 Mq (2.9 - 10 M©) from helium stars of 3 - 20 M© for 
/wR = 5 (/wR = 10; Table [2 cf. Fig. [HI). Our binary 
star models also give generally larger final masses for SN 
Ibc progenitors than those in WL99: Mf ~ 1.64 - 4.5 M© 
for /wR = 5 from Mzams =13 - 25 M©, compared to 
Mf ~ 2.0 - 3.4 M© in WL9^. 

5.2. Helium 

Figure [To] shows the amount of helium (Afnc) as a func- 
tion of the final mass in the SN Ibc progenitor models 
with /wR = 5 at solar metallicity. In general, both 
binary and helium star models predict large amounts 
of helium in the envelope compared to the results of 
WLW95 and WL99 except for relatively low mass sys- 
tems of Mne.i ^ 3.0 M©. Note that only rather massive 
progenitors with Mf > 5.5 M© can be significantly he- 
lium deficient (i.e., A/ho < 0.5 M©). 

For the systems with Mf < 3.0 M©, the role of Case 
ABB/BB mass transfer becomes important. In partic- 
ular, the small amount of helium (Afne ^ 0.18 M©) in 
Seq. 3 shows that helium deficiency can also be achieved 
by the so-called Case BB/ABB mass transfer from rela- 
tively low mass helium cores (AfHe,i < 3 M©). Therefore, 
we expect two different classes for Type Ic progenitors at 
solar metallicity: one with Mf < 2.0 M© and the other 
with Mf > 5 5 M© (see also iWellstein fc Langerl 119991: 
IPols fc Dewi1 [2002l. 

Obviously, the upper limit of Mf for potential Type Ic 
progenitors increases with a smaller mass loss rate: with 
/wR = 10, we have Mf > 8.9 M©. 

® Here, the lower end of Mf is determined by the detailed history 
of Case ABB/BB transfer, which depends on uncertain parameters 
such as the semi-convection efficiency. 



5.3. Hydrogen 

In our binary star models. Case AB or Case B mass 
transfer does not completely remove hydrogen from the 
primary stars, as shown in Fig. [Jj The resulting WR 
mass loss is therefore weaker than that from the corre- 
sponding pure helium star. This explains the fact that 
for a given WR mass loss rate, the binary star models 
with Mf > 3.0 M© predict somewhat larger A/hc than 
the single helium star models do as shown in Fig. [TOl 

More importantly, the binary star models show that, 
even aX Z k, Zq , small amounts of hydrogen can be re- 
tained up to the pre-supernova stage for a certain range 
of the final mass. Fig. [TT] (see also Table [T] and Fig. [IJ 
shows the total mass of hydrogen in the primary stars 
at neon/oxygen burning phase. The time span from this 
stage to core collapse is supposed to be less than about 10 
yr. The WR winds mass loss rate from the primary star 
models shown in the figure is about 10~^ M© yr~^. Some 
of the primary stars in this sample are still undergoing 
Case ABB or BB mass transfer (Seqs. 5, 8, 9, 13, & 21 
and many of the SMC models; see Table[T]), but the mass 
transfer rate is only about 10~^ M© yr~^. Therefore, the 
amount of hydrogen at core collapse should remain close 
to the values given in the figure. 

Fig [m indicates that the presence of hydrogen is only 
expected for 3.0 < Aff[M©] < 3.7, at solar metallicity 
with /wR = 5. This results from two different reasons for 
different mass ranges. For Mf < 3.0 M©, primary stars 
expand during the carbon burning phase to much larger 
radii than more massive ones do. The resulting mass 
loss rates via Case ABB or BB transfer thus become large 
enough to completely remove hydrogen from the primary 
stars by the time of core collapse. For Mf > 3.7 M©, on 
the other hand, WR winds are rather strong and can 
remove hydrogen from the primary stars during the core 
helium burning phase. 

The upper limit of Mf for the presence of hydrogen thus 
depends on the adopted mass loss rate. As shown in the 
figure, with /wr = 10, it increases to about 4.5 M© at 
solar metallicity. At SMC metallicity, the stellar wind 
effect is less significant and even rather massive super- 
nova progenitors of about 7 M© can retain fairly large 
amounts of hydrogen (~ 10~^ M©). We discuss implica- 
tions of this result for supernova types in Sect. 16.41 

5.4. Radius 

Relatively low mass helium stars usually experience a 
rapid expansion of the envelope during the core carbon 
burning phase. Interestingly, our new models predict 
much larger radii at the pre-supernova stages than what 
the WLW95 models do, as shown in Fig. [121 The first 
reason f or this difference is the u pdated OPAL opacity 
table bv llglesias &: Rogers I ()1996[ ). which causes a strong 
iron bump at around logT = 5.3. Another reason is the 
smaller WR mass loss rate adopted in the present study, 
and the resultant thicker helium envelope in our mod- 
els. The radius is further affected by the complication 
of binary interaction: the presence of a thin hydrogen 
layer in some of the binary models results in a more ex- 
tended envelope than in the corresponding single helium 
star models. It is also worth noting that the radii of the 
binary star models at Mf < 3 M© are smaller than those 
of the single helium star models. This is because the pri- 
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Fig. 12. — Tlie predicted radii of SN Ibc progenitors as a func- 
tion of the final mass. Tiie filled circles and triangles denote the 
results of our binary star models at Z = Zq for /wR = 10 and 5, 
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from mass-losing single helium star models at Z = Zq are marlsed 
by the dashed {/wR = 10), dotted (/wR = 5) and dashed-three- 
dotted (WLW95) lines. The thin dashed-dotted line gives the solar 
radius. 

mary stars in these sequences are filling the Roche-lobe, 
and the envelope cannot expand beyond it. 

The metallicity effect is somewhat subtle. The hy- 
drogen and helium layers in SN Ibc progenitors become 
thicker due to the reduced mass loss rate for lower metal- 
licity, while the opacity due to metals becomes smaller. 
As these two effects compensate each other, the radii of 
our SN Ibc progenitor models at both SMC and solar 
mctallicities are found to be similar. Finally, it should 
be noted that the radius depends on the final mass. In 
general, a less massive progenitor tends to have a larger 
radius at the presupernova stage. This might make mix- 
ing of nickel into the helium rich layer induced by the 
Ray leigh- Taylor instability during the supernova explo- 
sion more efficient for a less massiv e SN Ibc progenitor , 
as implied by the simulations of H achisu et al. I (|1991[ ) 
and lJoggerst et aTT (2009^ resulting in more prominent 
helium lines in the spectra (see Sect. [B] below). 

5.5. Winds 

Fig. [13] shows the mass loss rate due to WR winds 
at the central neon/oxygen burning phase from different 
SN Ibc progenitor models. The mass loss rate ranges 
from 10~^ Moyr^^ to ^ 10^^ M0yr~^, depending on 
the final mass, at solar metallicity. The corresponding 
wind velocity (i.e., the escape velocity) changes from ~ 
200 km s-^ to ~ 2400 km a'^. 

The circumstellar interaction of such winds from pri- 
mary stars in binary systems are supposed to be very 
complex because of the orbital motion and the wind- 
wind interaction with the secondary star. Some of the 
primary star models shown in the figure are still under- 
going Case ABB or BB mass transfer, which may add 
to the complexity. Therefore, the nature of the circum- 
stellar materials around SN Ibc progenitors in a binary 

''' However, it should be kept in mind that the degree of mixing 
due to the Rayleigh-Tayler instability may depend not only on 
the stellar s tructure, but also on explosion en ergy and directional 
asymmetry iHammer. Janlta fc Mueller 1120101) 
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system may not follow the simple relation of p oc r^^ 
that is expected for single WR star progenitors. 

6. DISCUSSION 

We have presented new evolutionary models of mas- 
sive close binary stars, considering tidal interaction, and 
transport of angular momentum and chemical species 
due to rotationally induced hydrodynamic instabilities 
and the Spruit- Tayler dynamo. We have investigated 
the redistribution of angular momentum in the primary 
star. Although mass transfer and tidal interaction can 
significantly affect the evolution of the rotation velocity 
of the primary star on the main sequence, the amount of 
angular momentum retained in the core in the late evolu- 
tionary stages is rather insensitive to the previous history 
of such binary interactions because of the self-regulating 
nature of the Spruit-Tayler dynamo. 

We have also calculated non-rotating, mass-losing sin- 
gle helium star models and compared them with our pri- 
mary star models in binary systems. Our models adopt 
a much lower WR mass loss rate than in the previous 
studies by WLW95 and WL99, and predict some new 
important properties of SN Ibc progenitors accordingly. 
The following discussions are based on the models with 
/wR = 5, unless otherwise specified. 

1. The final masses of SN Ibc progenitors in bi- 
nary systems at Z ~ Zq are not limited to 
1.5 Mq < Mf < 4 Mo as predicted by WLW95 and 
WL99, but a more wide range Mf is expected (i.e., 
1.5 Mq < Mf < 7.1 Mq from Mi„it ~ 12. ..60 M©; 
see Fig. [T4|l. 

2. At Z ~ Zq, significant deficiency of helium (A/ho < 
0.5 Mq) is expected for Mj/Mq > 5.5 and for 1.5 < 
Mf/MQ < 2.0 (Fig. [ini). Rather large amounts of 
helium up to 1.5 Mq are expected for the other 
final mass range (i.e., 2.0 < M{/Mq < 5.5 Mq). At 
Z ~ Zsmc , no such helium deficient SN progenitors 
are expected for the considered initial masses (16 — 

40 Mq). 

3. A thin layer of hydrogen with Mh — 10^* — 
10~^ Mq is predicted for SN Ibc progenitors with 
3.0 Mq < Mf < 3.7 Mq at Z ~Zq, and 3.0 Mq < 
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Mf < 8 Mq at Z 
Table [J). 



^SMC: 



respectively (Fig. (TTJ 



4. Most SN Ibc progenitors with Aff < 5 Mq rapidly 
expand during core carbon burning, resulting in 
i? =~ 4.0 — 30 Rq at the presuprnova stage. This 
is much larger than found in WLW95 (Fig. [T^ . 
Compact progenitors of i? < R© are only expected 
for a relatively high mass (Mf > 5.5 Mq at Z ~ Zq 
and Aff > 10 Mq at Z ~ Zsmc ; Fig.E]). 

The above results raise several important issues regard- 
ing observational consequences, as discussed below. 

6.1. Implications for energetic explosions powered by 
rapid rotation 

Our binary star models show that the mass transfer 
during helium core contraction (Case AB or Case B) in 
a close massive binary system cannot remove the hydro- 
gen envelope promptly enough to avoid the core braking 
due to the Spruit- Tayler dynamo during the helium core 
contraction phase. Compariso n of our binary star models 
with t he sing le star models bvlHeger. ^yooslev fc Spruiti 
(|2005D and lYoon. Langer fc Norman I (|2006D indicate 
that the amount of angular momentum retained in the 
core of the primary star at the presupernova stage should 
not be much different from those found in single star 
models if the Spruit-Tayler dynamo is adopted. I.e., a 
specific angular momentum of a few 10^^ cm^ s^^ in 
the innermost ^1.4 Mq at the presupernova stage is 
expected in both single and binary stars. This value 
is smaller by one or two orders of magnitude than 
what is necessary to make a long gamma-ray bursts by 
magnetar or coUapsar formation, or very energetic su- 
pernovae (hypernovae) powe red by rapid rotation and 
strong magnetic fields fe.g. iBurrows et al~ll2007[ ). al- 
though it may suffice to produ ce millisecond pulsars 
(iHeger. \yooslev fc SoruTtl 120051) . Together with the 
work bv iPetrovic et al. I (|2005al ). our results thus in- 
dicate that binary interactions with Case AB/B mass 
transfers at Z « Zq may not particularly enhance the 
production of strongly rotation-powered events like long 
GRBs or hypernovae. This is consistent with the obser- 
vational evidence that such events are rare compared to 
norinal core collapse events (e.g., iPodsiadlowski et al. I 
I200I iGuetta fc Delia Valiel 120071 ) . This also confirms 
the theoretical consensus that other evolutionary paths 
are needed to produce long GRBs associated with SN 
Ibc, such as the quasi-chemically homoge neou s evo- 
lution of a metal poor star dYoon fc Langer I 12005 



Yoon. Langer fc Norinan II2006I : iWoosley fc Heger I 



2006 



Cantiello et al. Il2007f ). tidal spin-up of a WR star in a 
very close binary s ystem with a neutron star or bla c k hole 
companion (e.g. I lzzard. Rami rez-Ruiz fc TouFI 120041: 
Ivan den Heuvel fc Y oon 2007j3, or binary evolution 
with Case C mass transfer with some specific co nditions 
(jBrown et al. |[2000l : IPodsiadlowski et al. IIMol) . 

6.2. Progenitor size 



^ A recent study using detailed stellar evolution models by 
IDetmers et al. I l|200SI ). however, seriously questions this possibil- 
ity 
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Fig. 14. — The predicted final masses of the primary stars in 
massive close binaries that undergo Case B mass transfer, as a 
function of the zero-age main sequence (ZAMS) mass, based on 
our binary and helium star models with /wR = 5. The absicssa 
is given in log scale. The numbers right above the abscissa denote 
the initial masses of the helium stars that the primary stars of the 
corresponding ZAMS masses would produce. The expected final 
outcomes of the primary stars according to different ZAMS masses 
are given by the labels right below the top: white dwarf (WD), type 
Ic supernova (Ic), type lb supernova with a thin hydrogen layer (lb 
with H) or type lib supernova (lib), and type lb supernova without 
hydrogen (lb without H). Here we assumed Mho < 0.5 Mq for 
SNe Ic progenitors. Note that the each boundary would shift to a 
higher value of Af2A.MS fo'' close binary systems with Case A mass 
transfer. 




Fig. 15. — Same as in Fig. [H but with /wR = 10. 

Larger radii of our SN Ibc progenitor models than 
those previously found should have consequences in 
shock break-outs and bolometric light curves. For in- 
stance, a shock break-out from a larger envelope would 
be marked by a lower photosperic temperature . Detailed 
comparison of numerical calculations with observational 
data may thus give st rong constraints on SNe Ibc p rogen- 
itor properties (e.g.. iCalzavara fc Matznei~ll2004[) . Re- 
cent discovery of t he X- ray outburst with SN 2008D by 
ISoderberg et al. 1 (12008') indeed suggests the usefulness 
of s uch a study for th e probe of supernova progenitors 
/e.g.lSoderberg et al. | l20 08: Che valier fc Fransson II200I 
iXu. Zou fc Fan I l2008inModiaz 6^1 I I2009D . foT^which 
our new models would provide ideal input. We will ad- 
dress this issue in a forthcoming paper. 

6.3. Presence of helium and implications for SNe Ic 
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Although the weak signature or no evidence of hehum 
in SNe Ic spectra may indicate the deficiency of hehum 
in their progenitors, it is not weh known how much he- 
hum can be hidden in the supernova spectra. It may 
also depend on the d egree of mixin g of n ickel into helium 
rich layers ([Wooslev fc Eastman! 119951) . If we assume 
0.5 Mq, for instance, as the maximum amount of he- 
lium allowed for hiding helium lines in SN spectra, our 
models indicate that most SNe Ic progenitors at solar 
metallicity should belong to two distinct classes in terms 
of both ZAMS and final masses, as summarized in Fig. [Til 
for the binary systems that undergo Case B mass trans- 
feiB. The same conc lusion was also drawn by WL99 and 
iPols fc DewTI (|2002[ ). But Pols & Dewi considered differ- 
ent types of binary systems (see below), and the finding 
of " two mass classes" for the final masses was not obvi- 
ous in WL99 due to the very high WR rate adopted in 
their study, although it was clearly seen for the ZAMS 
masses. 

If we assume that primary stars of 12.5 < 
-M'zams/Mq < 13.5 produce low-mass-class SNe Ic via 
Case BB mass transfer, about 62% of the SNe Ic from 
Case B systems should belong to the high mass class and 
the rest (^38%) to the low mass class, at Z « Zq (see 
Fig. [H]) . It is important to note that the two classes are 
produced by different mechanisms. The high mass class 
of SNe Ic progenitors (i.e., Mzams > 33 Mq in FiglM)) 
is a consequence of WR winds mass loss, while the low 
mass class results from Case ABB/BB mass transfer as 
discussed in Sect. 15. ll The ZAMS mass range for the low 
mass class may not be much affected by metallicity, while 
it should be widened with increasing metallicity for the 
high mass class. This leads to the conclusion that the 
low and high mass classes would dominate at low and 
high metallicities respectively. It should also be noted 
that the final mass range for the low mass class may not 
change much for different metallicities, while it may de- 
crease with increasing metallicity for the high mass class, 
due to the increasing WR winds mass loss rates, as im- 
plied by the resuh of WL99. 

In the massive close binary systems considered in this 
paper, the primary star masses become much smaller 
than those of the secondary stars when Case ABB or 
Case BB mass transfer begins. If the companion star 
mass were lower than the helium star in a close binary 
system, the mass transfer rate should become higher than 
in the systems of the present study. For example, if a 
helium star is located in a very short period binary sys- 
tem (P ^ 1 day) with a less massive companion (e.g. a 
neutron star), mass transfer from the helium star may 
occur rapidly enough to make a helium-deficient carbon 
star, e v en for Mhp.; ~ 6.0 M f^ as shown by Pols & Dcwi^ 
(pool . iDewi et al. I (po5l ) and llvanova et al. I (|2003[ ). 
The final masses of such SN Ic progenitors may range 
from 1.5 Mf7) to 3.0 Mf?) . Thi s scenario was also sug- 
gested bv lNomoto et al. 1 ()1994D to explain the fast light 
curve of Type Ic SN 19941. The ZAMS mas s of such a SN 
Ic should be in the range of 12 ~ 20 M© (jPols k, Dewi I 
[2002) . However, such close helium star plus neutron star 

^ For the systems with Case A mass transfer, the each boundary 
should move to the right in the figure, but the parameter space 
explored with our model grid is not large enough to determine it 
quantitatively 



systems are supposed to rarely form, and might not con- 
tribute much to the population of SNe Ic, compared to 
the systems considered in this study. 

Our results should have several observational conse- 
quences. As mentioned above, the population of SNe 
Ic should be dominated by the high mass class for 
Z > Zq. They would have higher ZAMS and fi- 
nal masses than those of typical SNe lb progenitors 
fcf. iKellv. Kirshner fc Pahre 1 120081: [Anderson k, James I 
l2008i r Given that the parameter space for the high 
mass class SNe Ic may become larger with a higher 
WR mass loss rate, the number ratio of SNe Ic to SNe 
lb, and that of high mass class SNe Ic to low mass 
cla ss SNe Ic should increase with increasing metallicity 
(cf.lPri eto. Stanek fc Beacom '12 0081 : [Anderson fc James I 
120091 iBoissicr fc Prantzos 2009). The existence of the 
two mass classes of SNe Ic progenitors may be re- 
lated to some aspects of the observational diversity of 
SNe Ic. For example, SNe Ic of the low mass class 
is likely to be characterized by rathe r fast declining 
light c urves and low luminosities fcf. Hwamoto et al. I 
'1994'; Richardson. P.. Branch. P.. fc Baron. E. I 120061 
Young et al. 2009), implying that the observed popu- 
lation of SNe Ic is likely to be baised to high-mass-class 
SNe Ic. 

We should also note the huge difference of the binding 
energy between the two classes. The binding energy of 
the envelope above 1.4 M© in the SN Ic progenitor star 
model of Seq. 3 (Mf = 1.64 Mg) is only about 10'*9 erg, 
while it should be one or two orders of magnitude higher 
for a SN Ic progenitor with Mf > 5.5 Mq (see Table [2]). 
As a consequence, the energetics of SNe Ic might be sys- 
tematically different for the two different classes. 

On the other hand, the assumption of Mho < 0.5 M© 
for SN Ic progenitors leads to a ratio of type Ic to lb 
supernova rate (Ic/Ib ratio )of about 0.4 from binary 
systems at solar metallicity Ej- This appears in contra- 
diction with recent observations t hat indicate r ather a 
high Ic/Ib ratio of about 2.0 (e.g. iSmartt ll2009h . This 
discrepancy would become even larger with /wR = 10, as 
implied by Fig. 1151 This raises a question on the nature 
of SN Ic progenitors, and it should be kept in mind that 
we still do not fully understand what distinguishes SN Ic 
progenitors from those of SN lb. 

A recent work by Pcs sart et al. I (|2010( ) indicates that 
the mass fraction of helium in the outermost layers (1^), 
rather than the total mass of helium, may be more rele- 
vant for the presence of helium lines in supernova spectra. 
Specifically, it is shown that if helium is well mixed with 
CO material such that becomes less than about 0.5, 
helium lines are not seen in early time spectra, despite 
rather a large total amount of helium (Mhc — 1.0 M©), 
if non-thermal excitation is absent. In our progenitor 
models, such a small is realized only for Mf > 5.5 M© 
at solar metallicity (with /wr — 5). This is not dif- 

At SMC metallicity, the SN Ibc progenitor models of initial 
masses of 16 - 40 Mq have 3.9 < Mf < 12 Mq with Mho > 1-2 Mq 
(see Table (TJ . This implies that the low-mass class SNe Ic would 
predominantly occur in binary systems at this metallicity. If we 
assume stars with 12.5 < M^^it < 13.5 would produce low-mass 
class SNe Ic as in the case of solar metallicity, the SN Ic/SN lb 
ratio would be about 0.1 at SMC metallicity. The exploration 
of the exact mass range for the low-mass class SNe Ic is a time- 
consuming task, and we plan to investigate this in near future. 
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ferent from the above-discussed mass limit for having 
MB_e ^ 0.5 Mq, implying that the initial mass range 
for SN Ic progenitors would not change much even if we 
adopted 1^ as a criterion, at least for the high mass class. 
On the other hand, we have — 0.98 in the primary 
star of Seq. 3 at carbon exhaustion while the total mass 
of helium is less than 0.2 Mq. This imples that the ini- 
tial mass range for the low mass class SN Ic progenitors 
might be affected if the condi tion oi < 0.5 for S N Ic 
progenitors were applied. But 'Dess art et al. I (|2010 D did 
not yet calculate such a low mass SN progenitor model 
(A/f < 2.0 M0), and their analyses were limited to early 
times of supernovae. It remains to be an important sub- 
ject of future work to systematically investigate which 
types of supernova progenitors would lead to the pres- 
ence or absence of helium lines in the supernova spectra 
at different epochs, including the effect of non-thermal 
excitation. Therefore, the above discussion on Type Ic 
progenitors based on the total amount of helium should 
only be considered indicative at this stage. 



6.4. Presence of hydrogen 

It is interesting that, at Z w Zq, the presence of a 
thin hydrogen layer is only expected for a limited range 
of the initial/final mass of SN lb progenitors, as shown 
in Figs [TTJ and [141 The detection of hydrogen ab- 
sorption lines at high velocity has b een indee d reported 
m many SNe lb (e.g.. IDeng et al. I [2000: Bra nch et al. I 
[2OOI lElmhamdi et al. II2006D . in favor of our model pre- 
diction for the presence of a thin hydrogen layer in SNe lb 
progenitors. This might provide a strong constraint for 
the progenitor masses of observed SNe lb, in principle. 

Note also that explosions of such helium stars with 



thin hydyrogen layers could be recognized as SN lib 
rather than lb, if hydrogen lines were detected short 
after supernova ex plosion, e.g ., as in the case of SN 
2008ax ( Cliornock et al. 1120101) a nd as recently di s cusse d 
by iSpencer fc Baron I (|2010^ and IDessart et al~l (|2010D . 
The radii of these progenitor models range from ~ 
10^^ cm to 10^^ cm. They may corredpond to the 
"compact " category of SN lib progenito rs, which is dis- 
cussed in iChevalier fc: Soderberg I (|2010[ ). The relatively 
low ejecta masses of such SNe lib are consitent with our 
model predictions. 

On the other hand. Case C mass transfer can also 
leave helium cores covered with small amounts of hy- 
drogen envelope. As the life time of such stars made via 
Case C mass transfer should be rather short, they can 
retain much more hydrogen (Mh > 0.1 Mq), than what 
is predicted from our binary models with Case AB/B 
mass transfer. Such a star r nay eventually explode as a 
SN lib like SN 199 3J (e.g. IPodsiadlowski et al. I [l99l 
iMaund et al. I I2004D , with a much extended envelope 
(- 10^3 _ 

Therefore, the two categories of SN lib progenitors ac- 
cor ding to their sizes, which has b een recently suggested 
bv IChevalier fc Soderberg"! ([2010l ). may be understood 
within the framework of binary evolution; SNe lib of the 
compact type may be produced via Case AB/B mass 
transfer (especially at Z < Zq), and SNe lib of the ex- 
tended type via Case C mass transfer. 

This work is supported by the DOE SciDAC Program 
(DOE DE-FC02-06ER41438), the NSF grant (NSF- 
ARRA AST-0909129), and the NASA Theory Program 
(NNX09AK36G). We are grateful to Luc Dessart for use- 
ful discussions. 
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TABLE 1 

Properties of the computed sequences 



No. 


z 


Ml i 


-''^2.1 


Pj 


/wR 


Case 




iwx f 




Mil 




^ J1.4 ^ 


Fate 


1 


0.02 


12 


8 


3.0 


5 


B+BB 


57.9 


1.40*= 


1.21*= 


0.17= 


0.0 


0.33 


ONeMg WD 


2 


0.02 


12 


11 


4.0 


5 


B+BB 


104.4 


1.48 


1.24<= 


0.20= 


0.0 


0.35 


ONeMg WD 


3 


0.02 


13 


11 


5.0 


5 


B+BB 


123.3 


1.64<= 


1.43*= 


0.18= 


0.0 


0.22 


SN Ic 


4 


0.02 


14 


12 


3.0 


5 


A+AB+BB 


118.5 


1.33= 


1.09= 


0.22= 


0.0 


- 


ONeMg WD 


5 


0.02 


14 


12 


5.0 


5 


B+BB 


30.7 


2.97 


1.66 


1.24 


1.9(-4) 


0.25 


SNIb 


6 


0.02 


16 


14 


2.0 


5 


A: Contact 














7 


0.02 


16 


14 


3.0 


5 


A+AB+ABB 


101.8 


1.54<= 


1.33= 


0.17= 


0.0 


0.39 


ONeMg WD 


8 


0.02 


16 


14 


4.0 


5 


B+BB 


26.2 


3.66 


2.05 


1.47 


4.5(-3) 


0.24 


SNIb 


9 


0.02 


16 


14 


5.0 


5 


B+BB 


33.7 


3.65 


2.04 


1.47 


5.0(-3) 


0.25 


SNIb 


10 


0.02 


18 


12 


3.0 


5 


A+AB+ABB 


27.9 


2.66 


1.58 


1.01 


0.00 


0.26 


SNIb 


11 


0.02 


18 


12 


3.0 


10 


A+AB+ABB 


27.3 


2.74 


1.59 


1.08 


0.00 


0.26 


SNIb 


12 


0.02 


18 


12 


5.0 


10 


B: Contact 
















13 


0.02 


18 


17 


3.0 


10 


A+AB+ABB 


36.2 


3.03 


1.68 


1.27 


7.9(-4) 


0.25 


SN lb 


14 


0.02 


18 


17 


4.0 


10 


A+AB 


29.7 


3.79 


2.14 


1.49 


1.5(-3) 


0.25 


SN lb 


15^ 


0.02 


18 


17 


4.0 


10 


A+AB 


24.0 


3.97 


2.27 


1.53 


1.0(-3) 


0.26 


SNIb 


le'' 


0.02 


18 


17 


4.0 


10 


A+AB 


50.0 


3.80 


2.16 


1.50 


2.2(-3) 


0.06 


SNIb 


17" 


0.02 


18 


17 


4.0 


10 


A+AB 


25.2 


3.84 


2.18 


1.50 


1.4(-3) 


0.26 


SNIb 


18'! 


0.02 


18 


17 


4.0 


10 


A+AB 


30.6 


3.73 


2.14 


1.43 


0.00 


3.57 


SNIb 


19 


0.02 


18 


17 


5.0 


3 


B 


33.1 


3.73 


2.33 


1.23 


0.00 


0.25 


SNIb 


20 


0.02 


18 


17 


5.0 


5 


B 


32.4 


4.04 


2.45 


1.4 


0.00 


0.33 


SNIb 


21 


0.02 


18 


17 


5.0 


10 


B+BB 


31.5 


4.41 


2.51 


1.68 


9.9(-3) 


0.26 


SNIb 


22 


0.02 


18 


17 


6.0 


10 


B 


39.3 


4.39 


2.56 


1.62 


4.0(-3) 


0.26 


SNIb 


23 


0.02 


25 


19 


6.0 


10 


B: Contact 














24 


0.02 


25 


24 


2.0 


10 


A: Contact 
















25 


0.02 


25 


24 


3.0 


3 


A+AB 


22.7 


3.70 


2.46 


0.98 


0.0 


0.24 


SNIb 


26 


0.02 


25 


24 


3.0 


5 


A+AB 


22.4 


4.33 


2.80 


1.30 


0.0 


0.25 


SNIb 


27 


0.02 


25 


24 


3.0 


10 


A+AB 


21.3 


5.07 


3.17 


1.67 


0.0 


0.25 


SNIb 


28^^ 


0.02 


25 


24 


3.0 


10 


A+AB 


18.9 


5.08 


3.19 


1.66 


0.0 


0.32 


SNIb 


29 


0.02 


25 


24 


4.0 


5 


A + AB 


21.5 


4.45 


2.91 


1.22 


0.0 


0.26 


SNIb 


30 


0.02 


25 


24 


6.0 


10 


B 


27.4 


6.49 


4.45 


1.63 


0.0 


0.39 


SNIb 


31 


0.02 


60 


40 


7.0 


3 


A 


16.8 


4.95 


3.70 


0.25 





0.24 


SNIc 


32 


0.004 


16 


12 


3.0 


5 


B+BB 


64.75 


3.91 


2.22 


1.54 


1.6(-2) 


0.24 


SNIb 


33 


0.004 


16 


14 


3.0 


5 


B+BB 


19.6 


3.90 


2.21 


1.53 


1.6(-2) 


0.24 


SNIb 


34 


0.004 


16 


14 


5.0 


5 


B+BB 


21.8 


3.84 


2.19 


1.51 


1.2(-2) 


0.24 


SNIb 


35 


0.004 


18 


12 


5.0 


5 


B+BB 


14.4 


4.64 


2.76 


1.68 


1.7(-2) 


0.31 


SNIb 


36 


0.004 


18 


12 


8.0 


5 


B+BB 


24.33 


4.56 


2.68 


1.67 


1.5(-2) 


0.33 


SNIb 


37 


0.004 


18 


17 


3.0 


5 


A+AB 


18.4 


4.42 


2.55 


1.67 


2.6(-2) 


0.26 


SNIb 


38 


0.004 


18 


17 


3.0 


10 


A+AB 


14.4 


4.58 


2.61 


1.77 


2.7(-2) 


0.26 


SNIb 


39 


0.004 


18 


17 


6.0 


5 


B+BB 


27.0 


4.57 


2.71 


1.65 


1.6(-2) 


0.27 


SNIb 


40 


0.004 


25 


12 


3.0 


5 


A; Contact 














41 


0.004 


25 


12 


6.0 


5 


B: Contact 
















42 


0.004 


25 


19 


3.0 


5 


A+AB 


10.2 


7.09 


4.87 


2.03 


6.5(-3) 


0.32 


SNIb 


43 


0.004 


25 


24 


3.0 


5 


A+AB 


13.0 


7.31 


5.05 


2.07 


8.7(-3) 


0.28 


SNIb 


44 


0.004 


25 


24 


6.0 


5 


B: Contact 














45 


0.004 


40 


30 


4.0 


5 


A+AB 


9.31 


12.0 


9.42 


1.24 


0.0 


0.56 


BH 



Each column has the following meaning. Z: the adopted metallicity, Mi^\: the initial mass of the primary component in units of Mq, M^^i- the 
initial mass of the secondary component in units of Mq, P\: the initial orbital period in units of day, /wr: the correction factor for the WR winds 
mass rate given in Eq. Case: the mass transfer ease, Pp: the orbital period in units of day, at the end of calculation (mostly at neon burning 
in the primary star), Mi f : the mass of the primary star at the end of calculation in units of Mq, Afco,f: the CO core mass of the primary star 
at the end of calculation in units of Mq, Mhc: the helium mass in the envelope of the primary star at the end of calculation in units of Mq, Mr: 
the hydrogen mass in the envelope of the primary star at the end of calculation in units of Mq, < ji,4 >: the specific angular momentum in the 
innermost 1.4 Mq of the primary star at the end of calculation in units of lO^^cms"^, Fate: the expected final fate of the primary star. 

^ /sync - 10' 
^ /sync - 0.01 

Taync accordiug to Zahn 
Non-magnetic model 

^ The values are lower limits for Mi^f and Mco,f j 3.iid upper limits for i\/Hc since the calculation was stop long before the neon burning phase 
while the Case BB/ABB mass transfer phase was not finished. 
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TABLE 2 

Properties of the computed single helium stars at 
Z = 0.02 





Mhe i 


Mf 


-'''^ CO, I 




R 


-R95 


Et 


5 


2.8 


2.73 


1.51 


1.16 


163.54 


5.08 


0.01 


5 


3.0 


2.91 


1.60 


1.22 


33.40 


1.79 


0.03 


5 


4.0 


3.72 


2.15 


1.43 


8.01 


0.73 


0.18 


5 


6.0 


4.23 


2.77 


1.19 


4.19 


0.38 


0.25 


5 


8.0 


4.95 


3.40 


0.85 


3.07 


0.30 


0.29 


5 


10.0 


5.49 


3.93 


0.54 


1.59 


0.22 


0.47 


5 


15.0 


6.44 


4.85 


0.31 


0.59 


0.18 


0.55 


5 


20.0 


7.09 


5.44 


0.22 


0.48 


0.16 


0.66 


10 


3.0 


2.95 


1.60 


1.27 


34.96 


1.84 


0.03 


10 


4.0 


3.84 


2.20 


1.49 


6.94 


0.70 


0.15 


10 


6.0 


4.93 


3.14 


1.55 


4.29 


0.51 


0.19 


10 


8.0 


6.06 


4.13 


1.41 


2.55 


0.35 


0.33 


10 


10.0 


7.01 


4.96 


1.12 


1.94 


0.25 


0.55 


10 


15.0 


8.89 


6.80 


0.47 


0.54 


0.20 


0.82 


10 


20.0 


10.33 


8.21 


0.31 


0.40 


0.16 


1.17 



Each column has the following meaning. fwK- the correction 
factor for the WR winds mass rate given in Eq. (Q, MHe,i' the 
initial mass of the helium star in units of Mq, Mf. the final mass at 
the end of calculation (i.e., at neon burning) in units of Mq , Mco,£ ■ 
the CO core mass at the end of calculation in units of Mq, Mue- 
the amount of Helium in the envelope at the end of calculation in 
units of Mq, R: the radius at the end of calculation in units of Rq , 
Rg^: the radius that encompasses 95 % of the total mass at the 
end of calculation in units of Rq, -Et,: the binding energy of the 
envelope above 1.5 Mq at the end of calculation in units of lO^^erg 



